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Abstract. We study electron densities of eigenfunctions of ato- 
mic Schrodinger operators. We prove the existence of p"'(0), 
the third derivative of the spherically averaged atomic density 
p at the nucleus. For eigenfunctions with corresponding eigen- 
value below the essential spectrum we obtain the bound p'"(0) < 
— (7/12)Z 3 p(0), where Z denotes the nuclear charge. This bound 
is optimal. 

o 
o 

1. Introduction and results 

q ■ In a recent paper |5. the present authors (together with T. Hoff- 

mann-Ostenhof (THO)) proved that electron densities of atomic and 
molecular eigenfunctions are real analytic away from the positions of 
< JL | the nuclei. Concerning questions of regularity of p it therefore re- 

mains to study the behaviour of p in the vicinity of the nuclei. A 
general (optimal) structure-result was obtained recently j2]. For more 
detailed information, two possible approaches are to study limits when 
approaching a nucleus under a fixed angle lj G S 2 , as was done in 
^ ■ [2], and to study the spherical average of p (here denoted p), which is 

mostly interesting for atoms. The existence of p'(0), the first derivative 
of p at the nucleus, and the identity p'(0) = — Zp(0) (see (jl.l2|) be- 
low) follow immediately from Kato's classical result ^2] on the 'Cusp 
Condition' for the associated eigenfunction (see also [T^], flOj). Two of 
the present authors proved (with THO) the existence of p"(0), and, for 
densities corresponding to eigenvalues below the essential spectrum, a 
lower positive bound to p"(0) in terms of p(0) in jUJ. In the present 
paper we prove the existence of p'"(0) and derive a negative upper 
bound to it (see Theorem II .2|) . A key role in the proof is played by the 
a priori estimate on 2nd order derivatives of eigenfunctions obtained 
in [H] (see also Remark 11.31 and Appendix iBl below) . Furthermore our 
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investigations lead to an improvement of the lower bound to p"(0) (see 
Corollary II. 7j) . The bounds on p"(0) and p"'(0) in terms of p(0) are 
optimal (see Remarks II .51 and \1.8\i . 

We turn to the precise description of the problem. We consider a 
non-relativistic iV-electron atom with a nucleus of charge Z fixed at 
the origin in IR 3 . The Hamiltonian describing the system is given by 

N y 1 

H = H N (Z) = J2{-^-r^)+ E u^T (!•!) 

j=l 1 Jl l<i<j<N 1 Jl 

The positions of the iV electrons are denoted by x = (xi, x 2 , ■ ■ ■ , scjv) G 
R , where ) denotes the position of the j'th electron 

in IR 3 , and Aj denotes the Laplacian with respect to Xj. For shortness, 
we will sometimes write 

H = -A + V(x), (1.2) 
where A = Ylj=i Aj * s the 3iV-dimensional Laplacian, and 

N 

v(x) = vw,z(x) = x;-r-r+ E i 1 ( L3 ) 

j=i 1 Jl i<i<j<N 1 n 

is the complete (many-body) potential. With Vj the gradient with 
respect to Xj, V = (Vi, . . . , Vat) will denote the gradient with respect 
to x. 

It is a standard fact (see e.g. Kato [TT]) that H is selfadjoint with op- 
erator domain T>(H) = W 2,2 (M. 3N ) and quadratic form domain Q{H) = 
W 1,2 (R 3JV ). 

We consider eigenf unctions of if, i.e., solutions ip e L 2 (R 3N ) to 
the equation 

Hip = Ei/j, (1.4) 

with Bel. To simplify notation we assume from now on, without 
loss, that ip is real. Apart from the wave function ijj itself, the most im- 
portant quantity describing the state of the atom is the one-electron 
density p. It is defined by 

N N „ 

P(x) =EPi( x ) = E / Wz,%)| 2 d% , (1.5) 
where we use the notation 



Xj = (xi, . . . , Xj-i, x j+ i, x N ) (1.6) 
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and 

dxj = dx\ . . . dxj^id,Xj + i . . . dxN (1.7) 
and, by abuse of notation, identify (xi, . . . , Xj-i, x, Xj + i, . . . , xn) and 

(x,Xj). 

We assume throughout when studying p that E and ip in ()1.4|) are 
such that there exist constants Co, 7 > such that 

|V>(x)| < C e~ 7|x| for all x G R 3N . (1.8) 

The a priori estimate [01 Theorem 1.2] (see also [01 Remark 1.7]) and 
(jl.8j) imply the existence of constants Ci,ji > such that 

|V^(x)| < d e" 7l|x| for almost all xGl M (1.9) 

Remark 1.1. Since if) is continuous (see Kato [T2*]). (J1.8)) is only an 
assumption on the behaviour at infinity. For references on the ex- 
ponential decay of eigenf unctions, see e.g. Froese and Herbst [7] and 
Simon [Tl|. The proofs of our results rely (if not indicated otherwise) 
on some kind of decay-rate for ip; exponential decay is not essential, 
but assumed for convenience. Note that (jl.8|) and (jl.9|) imply that 
p is Lipschitz continuous in M. 3 by Lebesgue's theorem on dominated 
convergence. 

In j2] we proved (with THO) that p is real analytic away from the po- 
sition of the nucleus. More generally, for a molecule with K fixed nuclei 
at . . . , R K , Rj G M 3 , it was proved that p G C w (R 3 \{i2i, . . . , R K })\ 
see also [3] and 0] . Note that the proof of analyticity does not require 
any decay of p (apart from ip G W 2 ' 2 (M. 3N )). That p itself is not ana- 
lytic at the positions of the nuclei is already clear for the groundstate 
of 'Hydrogenic atoms' (N = 1); in this case, Vp is not even continuous 
at x = 0. However, as was proved in [2], e z ^p G C 1,:L (1R 3 ). 

To obtain more information about the behaviour of the density at 
the positions of the nuclei one therefore has to study the regularity of 
other quantities, derived from p. 

One possibility is to study the function r i— > p(r,u) := p(ru) for 
fixed u) = A G S 2 (r = \x\); results in this direction were derived by 
the authors (with THO) in j2J Theorem 1.5]. In particular, for the case 
of atoms it was proved that for all w 6 § 2 , 

p(-,u) G C 2 ' Q ([0,oo)) for all a G (0, 1) , (1.10) 

and the 1st and 2nd radial derivatives were investigated at r = 0. 
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The main quantity studied in this paper is the spherical average of 

P, 

p(r) — / p{ruj) duj , r G [0, oo). (1-H) 

It follows from the analyticity of p mentioned above that also p G 
C^((0,oo)), and from the Lipschitz continuity of p in IR 3 that p is 
Lipschitz continuous in [0,oo). 

The existence of p'(0), the continuity of p ' at r = 0, and the Cusp 
Condition 

p'(0) = -Zp(0), (1-12) 

follows from a similar result for ^ itself by Kato |12| : see [IS], PD], and 
Remark 1.13]. 

To investigate properties of p and the derived quantities above it is 
essential that p satisfies a differential equation. Such an equation easily 
follows via (jl.4|) from 

- r 

V/ il>(x,± j )(H-E)il){x,1tj)dk j = 0. (1.13) 

This implies that p satisfies, in the distributional sense, the (inhomo- 
geneous one-particle Schrodinger) equation 

1 Z 

--Ap-—p + h = in R 3 . (1.14) 
The function h in ()1.14|) is given by 

N 
3=1 

f N f Z 

^R 3Ar " 3 i ^J^ j J9? N -^ \ x i\ 

N 

+ E / i ,hM*,x,)| 2 rfx, (1.16) 

^jJ*™-* \X-Xi\ 

+ /~2 I i ~ 1 IV'(^ ) %)| 2 ^ - Epj(x). 

i<k<e7^,k^e J ^ N -^ Xk ~ x ^ 

The equation ()1.14j) implies that the function p in p. 11)1 satisfies 

1 Z ~ 

__Ap p + /i = for rG(0,oo), (1.17) 

2 r 
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where A = + --f = 4-4 (r 2 -f ), and 

ar^ r ar r z ar ^ dr > ' 

h[r) = / h(ru)du. (1.18) 
is 2 

In [0J Theorem 1.11] information about the regularity of if) was used 
to prove that h G C°([0, oo)), and, using ()1.17j) . that consequently, 
p G C 2 ([0,cx))), and 

p"(0) = l{m + Z 2 p(0))- (1-19) 
Moreover, denote by a(H N (Z)) the spectrum of H N (Z), and define 

e := E%_ X (Z) - E , BP N _ X {Z) = inf a(H N ^(Z)) . (1.20) 
Then if e > Theorem 1.11], we have 

h(x) > ep{x) for all x G M 3 , (1.21) 
and so in this case, ()1.19j) implies that 

p"(0)>l(Z 2 + s)p(0)>lz*p(0)- (1-22) 

Our main result in this paper is the following. 

Theorem 1.2. Let ip G L 2 (IR 3Ar ) fee an atomic eigenf unction, H^(Z)ifj 
= Eip, satisfying ()1.8)) . u>zt/i associated spherically averaged density p 
defined by (JT3j) and ffTTTTTl . Let A fee denned fey (ITToT) - (fTToT) . (ITTKJ) . 

and /e£ £ fee gwen fe?/ fll.20|) . Let finally <pj(x, x,-) = e ~2' x 'ip(x, x,-) ; j = 
1,.,.,/V^ 

T/ien p G C 3 ([0, oo)), and 

p"'(0) = /J'(0)-|[M0) + Z 2 p(0)] (1.23) 

7 W r /• 

= - -Z 3 p(0)-^ZJ2 / iV^Cx.Opdx, (1.24) 

+ |(^(0, •), [tftf-l^ " 1) - ^(0, •)) i2(R 3--3 ) 

lfe>0, then 

p"'(0) < ~ ^{7Z 2 + 20e) p(0) < - ^Z 3 p(0) . (1.25) 

Remark 1.3. The existence of p^-*(0) for all k > 3 remains an open 
problem. The two main steps in the proof of Theorem II .21 are Proposi- 
tions and below. From the latter one sees that the existence of 
h^ k ~ 2 '(0) is necessary to prove existence of p^-'(O). In Proposition 11.61 
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the existence of h'(0) is proved and this result already heavily relies on 
the optimal regularity results for ip (involving an a priori estimate for 
second order partial derivatives of if)) obtained in [OJ (see also Appen- 
dix El below). 

Remark 1.4. Note that the inequality 

P'"(0) < - ^ 3 P(0) (1.26) 
follows from ()1.24j) as soon as ip is such that 

N 

V (V>(0, •), [H N -i(Z - 1) - E}^(0, •) W*- 3) > . (1.27) 
i=i 

There are cases where ()1.27|) holds even if the assumption e > does 
not. For instance when E is an embedded eigenvalue for the full op- 
erator H (e < 0), but non- imbedded for the operator restricted to 
a symmetry subspace. In particular, (|1.26|) holds for the fermionic 
ground state. 

Remark 1.5. Compare (jEOJ), (fT22|) . and Theorem Ol with the fact 
that for the ground state of 'Hydrogenic atoms' (N = 1), the corre- 
sponding density pi(r) = ce~ Zr satisfies 

pl k \0) = (-Z)%(0). (1.28) 

In fact, if (-A - Z/\x\)ijj n = E n ifj n , E n = -Z 2 /4n 2 , n e N, ijj n {x) = 
e~ ~2^(j) n (x), then (jl.24j) implies that the corresponding density p n sat- 
isfies 

= [" ^ + l ZE ]Pn{0) -4vrZ|V0 n (O)| 2 
12 L rr J 

For the ground state, i.e., for n — 1, Ex — — Z 2 /4, 0! = 1, this reduces 
to (IT^HD with k = 3. 

Furthermore, for s - states (zero angular momentum), we get that 
V0 n (O) = 0, since n is radial and C 1,a (see f|3.25|) below). Taking n 
large in ()1.29|) illustrates the quality of the bound (|1.25|) . 

The proof of Theorem 11.21 is based on the following result on h. Its 
proof is given in Section 13.11 

Proposition 1.6. Let ip e L 2 (M. 3N ) be an atomic eigenf unction, 
Hn(Z)iP = Eip, satisfying (jl.Sj) . and let h be as defined in ()1.15j) - 
(fTTToT) . LetueS 2 andh(r) = J &2 h{ru) duj . 
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Then both h and the function r \— > h(r,u) := h(ru) belong to 
C\[0,oo)). 

Furthermore, with (fj(x, x,) = e~^^ip(x, x,-), j — 1, . . . , N, 



Ho) 



h'(0) 



Z 2 N r r 

— p(0)+4tt^ / |V,^(0,x,)| 2 rfx, (1.30) 
+ (V>(0, •), [H N -i(Z - 1) - E]^{0, •) W"- 3) l , 

7 3 A N f 

- Zh(0) + — p(0) + ^zJ2\ |V^(0, Xi )| 2 ^ 



-^(O,.),^-!^-!)-^,-))^^) . (1.31) 

As a byproduct of (jl.30j) we get the following improvement of (|1.22|) . 

Corollary 1.7. Letip e L 2 (R 37V ) oe an atomic eigenf unction, H^(Z)ifj 
= Eip, satisfying (jl.Hjl . wtt associated spherically averaged density p 
defined by (|1.5jl and fll.lljl . Let e 6e given 6?/ (|1.2()jl . and assume e > 0. 
Taen 

P"(0)>|[^f + e]p(0)>|z 2 p(0). (1.32) 

Proof. Using the HVZ-theorem [T3*l Theorem XIII. 17], ()1.30|) provides 
an improvement of the bound (|1.21|) for r = to 

h(0) > + e]p(0) ■ (1-33) 

This, using (jOty . gives |HT2jl . □ 

Remark 1.8. For 'Hydrogenic atoms' (N = 1) (see Remark II. 5|) . 
(fTHH) and (OHjl imply 

3?(0) = ^ [5 + l]p„(0) + ^| V0„(O)| 2 , (1.34) 
6 n^ J 3 

which illustrates the quality of the bound ()1.32|) above (see also the 

discussion in Remark II. 5j) , and reduces to (jl.28|) with k = 2 for the 

ground state (n = 1, = 1). 

We outline the structure of the rest of the paper. In Section 
we use Proposition 11.61 and the equation for p (see (|1.17|) ) to prove 
Theorem 11.21 In Section |3] we then prove Proposition 11.61 This is done 
applying the characterization of the regularity of the eigenfunction ip 
up to order C 1,1 proved in [H] (see also Appendix |B] and Lemma 13. 9j) 
to the different terms in ()1.15|) - (jl.l6|) . 
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2. Proof of Theorem 11.21 

That p G C 3 ([0, oo)) and the formula (jl.23j) follow from Proposi- 
tion |23 below (with k = 1), using Proposition 11.61 and (|1.19j) . The 
formula (jl.24j) then follows from (jl.23j) and Proposition 11.61 

If e > 0, then the HVZ-theorem Theorem XIII. 17] implies that 

TV 

X>(0, •), [H N - X {Z - 1) - Ety(0, 0)^(^-3) (2.1) 

3=1 

N 

> £ y)<^(0, ■), "0(0, -))ia(R? w - 3 ) = » 

which, together with (jl.24j) . implies (|1.25jl . since p(0) = 47rp(0). 
It therefore remains to prove the following proposition. 

Proposition 2.1. Let if) G L 2 (M. 3N ) be an atomic eigenf unction, 
Hn(Z)i/j = Eif), satisfying (jl.HJ) . wift associated spherically averaged 
density p defined by (jl.5|) and (JHHJ, and Zet /i 6e as defined in (|1.15|) - 
dnnj) and (dUHJI- Let A; G N U {0}. 

If h & C fc ([0, oo)) t/ien p G C fc+2 ([0, oo)) ; and 

p(fc+2)( ) = -1_ [(A; + l)h {k \0) - Zp ( fe+1 )(0)] . (2.2) 

Proof. Let r > 0; multiplying ()1.17|) with r 2 , integrating over [5, r] for 
< 5 < r, and then taking the limit 5 j 0, using that h, p, and p' are all 
continuous on [0, oo) (see the introduction), it follows from Lebesgue's 
theorem on dominated convergence that 

p'(r) = A jT [ - Zp(s)s + h(s)s 2 } ds 

= -2Z f p(ra)ada + ^- [ h(s)s 2 ds . (2.3) 
Jo r Jo 

Using again Lebesgue's theorem on dominated convergence (in the form 
of Proposition 13.61 below) we get that for r > 0, 

p "(r) = 2[h(r) - [ [Zp'(ra) + 2h(ra)] a 2 da] . (2.4) 
Jo 

Since p G C 2 ([0, oo)) (see the introduction), ()2.4|) extends to r = by 
continuity and Lebesgue's theorem. This finishes the proof in the case 
k = 0. 

For k G N, applying Lebegue's theorem to ()2.4j) it is easy to prove 
by induction that if h G C fc ([0, oo)) then p G C h+2 ([0, oo)), and that 
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for r > 0, 

p( k+2 \r) = 2[h {k) (r)- f [Zp {k+1 \ra) + 2h {k \ra)]a k+2 da]. (2.5) 

Jo 

In particular, (|2.2|) holds. This finishes the proof of the proposition. □ 

This finishes the proof of Theorem 11.21 □ 
It remains to prove Proposition II .61 

3. Study of the function h 

3.1. Proof of Proposition [TTBl It clearly suffices to prove the state- 
ments in Proposition II .61 for each hj (j — 1, . . . , N) in ()1.16j) . Propo- 
sition 11.61 then follows for h by summation. We shall prove the state- 
ments in Proposition 11.61 for hi] the proof for the other hj is completely 
analogous. 

Recall (see (fjT5]) - fTTB]l ) that hi is defined by 

h\(x) = ti(x) — V\(x) + w\(x) — Epi(x), (3.1) 
h(x) = [ |V^,xi)| 2 dx 1; (3.2) 

- r z 

M%) = /2 rl^.xi)! 2 ^!, (3.3) 

^ J R 3N-3 \X k \ 

N l 
W l( X ) = / J l 1 r |^(x,Xi)| 2 rfXi 



k=2 



-3 \X — Xk\ 



+ V / Mx^dki, (3.4) 

Pi(x) = / \ip(x, xi)| 2 rfxi. (3.5) 



Here, x x = (x 2 , . . . , x^) G M 3Ar ~ 3 and rfxi = <ix 2 . . . dx^. 

We shall look at the different terms in (J3.2j) - ()3.5)) separately. The 
statements on regularity in Proposition 11.61 clearly follow from Propo- 
sition and Proposition 13. 31 below and the regularity properties of p, 
see the discussion in the introduction (in the vicinity of (jl.l(JJ) - (jl.l2|) ; 
see also Theorem 1.5] and [HI Theorem 1.11]). The formulae 1)1.30)) - 
()1.31j) follow from the formulae in Proposition 13.11 and Proposition 13.31 
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using ()3.1j) - (j3.5j) . and the fact that 

/ (|V il V'(0,x 1 )| 2 +(V, v _ 1 ^_ 1 (x 1 )- J B)|V'(0,x 1 )| 2 ) d± 1 

JR3N-3 L J 

= <^(0, •), [H N ^(Z - 1) - £#(0, •)) i2(K |f-3 ) • (3.6) 

For the definitions of Vn-i,z-i an d Hn-i(Z — 1), see (|1.3jl and 
respectively. 

Proposition 3.1. Let t>i and W\ be defined as in ()3.3|) - (|3.4|) and Zet 

w G S 2 . Define Vi(r) = L 2 Vi(ru) du, W\{r) = L 2 W\(ruS) du. 

Then Vi(- ,u),wi(- ,u),Vi, and Wi all belong to C 1 ([0, oo)) ; and 

v{(0)= -Z5i(0), (3.7) 
w{(Q) = -ZwtiO). (3.8) 

Remark 3.2. Similar statements hold for Vj and Wj (j = 2, . . . , N), 
and therefore, by summation, for v = J2f=i Vj, w = ^, =1 Wj. Compare 
this with ()1.12j) : p'(0) = —Zp(0); that is, three of the four terms in h 
(see (HISJ)) satisfy the 'Cusp Condition' /'(0) = -Z/(0). 

For the remaining term in (|3.1|) . we have the following. 

Proposition 3.3. Let t\ be defined as in (|3.2|) and let u 6 § 2 . Define 

Then both t\ and the function r i— > £i(r, cj) := ti(ru;) belong to 
C l ([0, oo)) . Furthermore, with (fi(x, Xi) = e^^'^ijXi), 

^i(O) = 4-Pi(°)+ 47r f / |Vi^i(0,X!)| 2 dxi (3.9) 

+ / iVx^^,^)! 2 ^ 

~ ~ Z 3 Air r /■ 

t/(0) = -Zt 1 (0) + — Pi(0) + — Z / iViPi^xOfdx! (3.10) 
{ |Vx^(0, Xl )| 2 + (^v-i,z-i(xi) - £)|^(0, Xl )| 2 ) dkx] , 

R 3iV-3 L J J 

with Vn~i,z-i given by (jl.3j) . 

Remark 3.4. Note that in the case of Hydrogen (N = 1), with 4>(x) = 

t'(0)+Zt(0) = |([^ + E]p H (0) + 4vr|V0(O)| 2 ) , (3.11) 
with V0 = for the ground state. 
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3.2. Integrals and limits. In the proofs of Proposition ^. H and Propo- 
sition 13.31 we shall restrict ourselves to proving the existence of the 
limits of the mentioned derivatives as r j (e.g., lim^o v[(r, u)). The 
existence of the limits of the difference quotients (here, v[(0,u) : = 
lim r | Vl ( r '")~ Vl (°) ) ; and their equality with the limits of the derivatives 
(i.e., lim r j_o v[(r, u>)), is a consequence of the following lemma, which is 
an easy consequence of the mean value theorem. 

Lemma 3.5. Let f : [a, b] -»• R satisfy f G C°([a,b]), f G C°((a,b)), 
and that lim e | f\ a + e) exists. 

Then lim e j_ ^ a+e ^~ :l '^ exists and equals lim e j f'(a + e). 

Verifying the two first assumptions in Lemma 13.51 in the case of the 
functions in Proposition ^. H and Proposition ^. 3l follow exactly the same 
ideas as proving the existence of the limits of the mentioned derivatives 
as r | (which follows below), and so we will omit the details. 

When proving below the existence of the limits of the mentioned 
derivatives as r j 0, we shall need to interchange first the differentiation 
j-, then the limit lim r | , with the integration / K 3jv_ 3 • • ■ dki (or, for the 
spherical averages, with J g2 J R 3jv- 3 • • • dk\duS). We shall use Lebesgue's 
Theorem of Dominated Convergence in both cases; in the case of 4- 
in the form of Proposition 13.61 below, which is a standard result in 
integration theory (see e.g. PQ). The dominant (in Proposition 13.61 the 
function g) will be the same in both cases. 

Proposition 3.6. Let I cR be an interval, A any subset o/R 771 , and 
f a function defined on Ax I , satisfying the three following hypothesis: 

i) For all AG/, the function x \— > f{x, A) is integrable on A. 

ii) The partial derivative df/d\(x, A) exists at all points in Ax I. 

iii) There exists a non-negative function g, integrable on A, such 
that \df/d\(x, A) | < g\x) for all (x, A) eAxI. 

Then the function F defined by 

F(X) = [ f(x,X)dx 

J A 

is differentiable in I , and 

F'(X) = J ^(x,X)dx. 

Remark 3.7. Note that one can remove a set B of measure zero from 
the domain of integration A, without changing the two integrals above; 
it is therefore enough to check the three hypothesis on this new domain, 
A' = A\B. Note that the set B of measure zero must be independent 
of A. 
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Note that in ()3.2|) - (j3.5j) we can, for x = tuj ^ 0, restrict integration 
to the set 

Si(u) = {x! G R 3N ~ 3 x k ^0, ^r^tu, x t ^x k , (3.12) 

L \Xk\ 

for k, t G {2, . . . , N} with k ^ £ } , 

since IR 3JV_3 \ S\(u) has measure zero. The set Si(u) will be our A'. 
From the definition of Si(ui) it follows that for any r > 0, uj G § 2 , and 
xi G <Si(w), there exists an e > and a neighbourhood C/ C Si(u) C 
M 3V- 3 of ii guch that y in is C°° on £ 3 (ru;,e) x U C M 3Ar . It 

follows from elliptic regularity [H] that ip G C°°(Bs(ru, e) x [/). In 
particular, if G : M 3iV — > R is any of the integrands in (J3.2)) - (|3.5|) . then, 
for all uj G § 2 , the partial derivative dG w /dr(r, Xi) of the function 
(r, Xi) i— > Xi) = G ll) (r , 1 Xi) exists, and satisfies 

BC 

—^(r,x 1 )=oj-\(V 1 G)(ruj,± 1 )] for all r > , Xl G SAu). (3.13) 
ar L J 

This assures that the hypothesis ii) in Proposition 13.61 will be satisfied 
in all cases below where we apply the proposition (the hypothesis i) is 
trivially satisfied). 

We illustrate how to apply this. Let 



9w{r) = g(r,uj) = / G(ru,-kx) dx x (= / G{ru),-k x ) dkx) . 

Interchanging integration and differentiation as discussed above, we 
have, for uj G § 2 fixed and r > 0, that 

gj(r) = [ to - [(V l G)(ruj,x 1 )] d&i. (3.14) 

To justify this, and to prove the existence of lim r ^o gJ{ r ) we will need 
to prove two things. 

First, we need to find a dominant to the integrand in (J3.14|) . that is, 
a function $ w G L^M 3 ^ 3 ) such that, for some Rq > 0, 



UJ 



(ViG)(rw,xi)) < ^(xi) for all r G (0,R ), x x G (3.15) 



This will, by Proposition 13.61 also justify the above interchanging of 4- 
and the integral ( (|3.15p ensures that the hypothesis iii) is satisfied). We 
note that, with one exception, whenever we apply this, in fact $ = <&,,, 
will be independent of uj G S> 2 , and therefore 



p3V-3 



X 
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Secondly, we need to prove, for all w £ § 2 and Xi G Si(u) fixed, the 
existence of 



lim 

r|0 



(ViG)(rw,X!) . (3.16) 



This, by Lebesgue's Theorem of Dominated Convergence, will prove 
the existence of lim r j_ gj(r). 
To study 

g{r) — I g{r ) uj)duj= i / G(ra>,Xi) dk\ duj (3-17) 

i§ 2 i§ 2 JIR3JV-3 

note that, by the above, the set 

P- 3 xS 2 )\(U [5i(o;)x{ W }f 
wes 2 

has measure zero, and that, also by the above, the partial derivative 
dG/dr(r, (xi,u;)) of the function 

(r, (xi, u)) i ► G(r, (x 1; u)) = G(ru, Xi) 

exists for all r > and (xl, a;) G U weS 2 [Si (a;) x {c<j}] . As noted above, 
the dominants $ we will exhibit below when studying g(r,cu) will (ex- 
cept in one case) be independent of a; G S 2 , and so can also be used to 
apply both Lebegue's Theorem on Dominated Convergence, Proposi- 
tion H03 and Fubini's Theorem on the integral in (J3.17)) . This implies 
that 

limo'(r)= / / (lima;- \(ViG)(ru,±i)] \ du dki, (3.18) 

as soon as we have proved the pointwise convergence of the integrand 
for all u G § 2 and Xi G «Si(a>), and provided the mentioned dom- 
inant. (In the last, exceptional case, we will provide a dominant 
$ G L 1 ^" 3 x § 2 )). 

In the sequel, we shall use all of this without further mentioning, 
apart from proving the existence of a dominant, and the existence of 
the pointwise limits of the integrands on Si(uS). Also, for notational 
convenience, we shall allow ourselves to write all integrals over IR 37V-3 
instead of over Si{uj). 

3.3. Additional partial regularity. For the existence of pointwise 
limits the following lemma will be essential; it gives detailed informa- 
tion about the structure of the eigenfunction if) in the vicinity of the 
two-particle singularity x = 0. The lemma is reminiscent of more de- 
tailed results obtained earlier, see Proposition 2], t3J Lemma 2.2], 
and Lemma 3.1]. 



14 S. FOURNAIS, M. HOFFMANN-OSTENHOF, AND T. 0. S0RENSEN 

We need to recall the definition of Holder-continuity. 

Definition 3.8. Let O be a domain in W 1 , I e N, and a G (0,1]. 
We say that a function u belongs to C k,a (Q) whenever u G C k (Q), 
and for all j3 G N n with \/3\ = k, and all open balls B n (xo,r) with 
B n (xo, r) C Q, we have 

\D^u(x) -D p u(y)\ 
sup 1 _ - yUn <C(x ,r). 

x,y£B n (x ,r),x^y \ x U\ 

When k = and a G (0, 1) we also write C Q (fi) = C°' Q (fi). 

Lemma 3.9. Let uj G § 2 and x^ 1 G iSi(co>). T/ien t/iere exists an open 
neighbourhood U C <Si(w) C R 3iV_3 of ~k\ and e > swc/i that 

?/>(x, xi) = e~ 2" |1 Vi( :r ) xi) with (3.19) 

d^ i ip 1 EC 1 ' a (B 3 (0,e) xU) for all a e (0,1), /3 G N 3A ^ 3 . (3.20) 

Proof. By the definition (J3.12j) of <Si(a;) there exists a neighbourhood 
[/ C C M 3Af ~ 3 of x° G and e > such that 

Xj- ^ x fc for j, k G {2, . . . , N} with j + k , (3.21) 

a?j ; ^ , 7^ x for j G {2, . . . , iV} 

for all (x, x 2 ,..., x N ) G 5 3 (0, e) x U C M 3JV . 

Make the Unsafe?/; = e~f |a Vi. Using l|L"3|» . ffTTjl . and that A x .(|x|) = 
2|x| _1 , we get that (fx satisfies the equation 

A(pi — Z-^-r ■ Vi^i + (^ + E — Vi)(fi = (3.22) 

with 

- Z - 1 1 

F 1 (x,x 1 ) = ^---+^ 1 -+ Yl I j" ' (3-23) 

j=2 1 jl fc=2 |X Xfe| 2<j<fc<iV ^ Xfc| 

where, due to ()3.21|) . 

V 1 eC co (B 3 (0,e)xU). (3.24) 

Since the coefficients in ()3.22|) are in L°°(I?3(0, e) x U)), this implies, 
by standard elliptic regularity [HI Theorem 8.36], that 

<pi G C^iBsfi, e) x U) for all a G (0, 1). (3.25) 
Let 7] G N 3Ar ~ 3 , 1 77 1 = 1; differentiating (|3~2^|) we get 

A^O-^-V^r) (3.26) 

= -[(f +^-Vi)(^ 1 )-(^ l y 1 )^ 1 ]. 
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By (f3~2~5jl and (jR~2Tty . the right side in belongs to C a (B 3 (0, e) x U) 

for all a G (0, 1) and so, by elliptic regularity, d^tpi G C 1,a (i?3(0, e)xU) 
for all a G (0, 1). An easy induction argument finally gives that 

S^pi G ^(^(O^) x [/) for all a G (0, 1), /3 G N 3N ~ 3 . 

□ 

3.4. Proof of Proposition [3~TT1 We will treat the two kinds of terms 
in (|3.4|) separately. 

For the first one, we make a change of variables. Assume without 
loss of generality that k = 2, and let y = X2 — ru, then 

1 



-3 \ruj — %2\ 



(ra;,xi) dxi 



(3.27) 



1 



R3iV-3 1 1/| 



(ro;,2/ + ru;,X3, 



XAf)| 2 dydxs . . . dx 



N 



For w G § 2 fixed, and r > 0, interchanging integration and differentia- 
tion as discussed above, we get, using (|3.27|) . that 



d 
dr 



1 



3 \ruj — X2 
2 



(ru;,xi)| 2 dxi 



tfj{u ■ (ViV> + V 2 V0} (ro;,j/ + ro;,a;3,...,a;jv)) (3.28) 



c/x 3 . . . cfe/v. 

Note that, using (jl.8j) . (jl.9|) . and equivalence of norms in M 37V , there 
exists positive constants C, ci, C2 such that 

2 ■ • (ViV> + V 2 V ; )}(™,1/ + ^,^3, • • -,xn) 



\y\ 



i 



(3.29) 



< Ce Cir e -c 2l(f> x 3i— j^jv)! 

Il/I 

which provides a dominant, independent of uj G S 2 , in the sense of 
(j3.15|) . uniformly for r G (0, i?o) f° r an y -Ro > 0. 

Writing ip = e~~^(^i, the integrand in (|3.28j) equals 

1 — r [ - Z|^| 2 + 2<fi e~ Zr (u ■ {Vi9?i + V 2 y?i})] (ra>, y + ru,x 3 , . . . , x N ), 

\y\ 

which has a limit as r | (for uj G S 2 and (y, a; 3 , . . . , xn) fixed) since 
(pi G C 1,OL by Lemma 13.91 This proves the existence of the limit of 
the integrand in (J3.28)) . and therefore of the limit as r [ of the 
derivative with respect to r of the first term in (|3.4|) . Note that since 
J s2 u (ko — 0, the terms proportional to w vanish by integration over 
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(y, £3, . . . , xn, to)- The limit of the derivative of the spherical average 
of this term is therefore (after setting x 2 = y) 

-Z I I |^(0, xi)| 2 rfici cfc^ . (3.30) 

is 2 \X 2 \ 

For the second term in ()3.4|) . assume without loss that k — 2, £ — 3, 
and write ip = e~ 2"' z 'y?i as before. Then, for u G § 2 fixed, and r > 0, 
interchanging integration and differentiation we get 



d 
dr 



-3 \x 2 - x 3 \ 
2 



(x, x x )| 2 diti 



ipi/ru, xi) |"c<j • Viipiru, Xi)l rfxi . (3.31) 

-3 \x 2 - x 3 \ 

1 [ - Z\^\ 2 + 2ip x e~ Zr (u ■ Vi^i)] (rw, x x ) dxi. 

-3 p2 — X3I 

(3.32) 

As above, (jl.8|) and (jl.9|) provides us with a dominant to the integrand 
in (|3.31|) in the sense of (|3.15|) . 

Also as before, the integrand in (J3.32)) has a limit as r j 0, since 
ipi G C ,a by Lemma 13.91 This proves the existence of the limit as 
r I of the derivative of the second term in (|3.4|) . Again, the term 
in ()3.32|) proportional to uo vanish when integrating over (xi,u;), since 
f s2 uidu) = 0. The limit of the derivative of the spherical average of 
this term is therefore 

-Z I I : |^(0,xi)| 2 dxida;. (3.33) 

J§2 JR3N-S \X 2 — X 3 \ 

This proves the existence of lim^o w{ (r, uj) (for any u G § 2 fixed), 
and of lim r j w{ (r). Furthermore, from (|3.3Uj) and (|3.33|) . 

N 

\imw{(r)=-z[ [ [V— 1 |V(0, Xl )| 2 dxida; 

r t° JR3N-3 \x k \i 

-zf f IV - r l ^(0, Xl )| 2 dx lC ?a; 

^ 2 L 2 <^ <JV Fk-a*N 



- Zwi(0). 
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The proof for v\ is similar, but simpler; we give it for completeness. 
For r > we get, by arguments as for w{, that 

N f 2Z 

v{(r,uj) = / - — | ip(ru>, xi) \uj ■ V\ip(ruj, xi)] dxi (3.34) 

^ JR3AT-3 \X k \ 

N r z 

= > / i— r [ - + 2<^i e" Zr (w ■ Vi^i)l (rw, x x ) dii. (3.35) 

One provides a dominant to the integrand in ()3.34|) in a similar way as 
for w\. We omit the details. Again, existence of the limit as r J, of 
the integrand in ()3.35|) is ensured by Lemma f3. 91 

The last term in (|3.35|) again vanishes when taking the limit r | 
and then integrating over (xj, u), since J g2 u duo = 0, and so 

N f f Z 
lim wi (r, cj) = — Zy \ \ - — - 1^(0, xi)| 2 dkidu 

= -Zvi(0). 

This finishes the proof of Proposition 13.11 □ 

3.5. Proof of Proposition We proceed as in the proof of Propo- 
sition Interchanging integration and differentiation as discussed in 
Section EIH we have, for w G S 2 fixed and r > 0, that 

t{(r,u)= ! u- [VidVVHCr^xOjdxi. (3.36) 

JR3JV-3 

Again, to justify this and to prove the existence of lim r iot{(r,u) 
we need to prove two things: The existence of the pointwise limit as 
r I of the integrand above, for to 6 S 2 and xi G fixed, and the 

existence of a dominant in the sense of ()3.15|) . 

Pointwise limits. We start with the pointwise limit. We allow our- 
selves to compute the integrals of the found limits, presuming the dom- 
inant found. We shall provide the necessary dominants afterwards. 

Recall from (the proof of) Lemma 1331 that <pi = e^^ip satisfies the 
equation 

A(pi _ ■ V m + (f + £ - V l )i Pl = (3.37) 
\x\ 

where 

V 1 eC°°(B 3 (p,e)xU) 
for some e > and U C Si(u>) C R 3Af ~ 3 some neighbourhood of xi. 
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Using Lemma ET31 we get that A^tpi G C x ^{B z {Q,e) x U) for all 
a G (0, 1). From this and ()3.37j) follows that for any Xi G Si(ou) fixed, 
and some e > 0, the function x h- > y>i(x, x x ) satisfies the equation 

A x(pi - Z^- ■ V^i = -A* (p! -(*+E- Vjcp! (3.38) 
\x\ 

= G kl , G kl G C x ' a (B 3 (0, e)),ae (0, 1) . (3.39) 

Furthermore, from Lemma f3 . 91 f with (5 = 0) it follows that Vi<^i(0, Xi) 
is well-defined. Note that for c G M 3 the function 

v(x) = -\x\ 2 (c ■ u>) = -\x\(x ■ c) (3.40) 

solves A x v = c ■ uj. Therefore, with c = ZVi<^i(0, xi), the function 
u = Mx x = (fx( ■ ,xi) — v satisfies the equation (in x, with Xi G Si(u) 
fixed) 

A x u(x) = Z— ■ (Vi<^i(x, x x ) - Vi^i(0, xi)) 

- Aa^ifoxO -(Zl + E- V x )<pi{x,x. x ) = g^(x). (3.41) 

By the above, and Lemma TA. II in Appendix 1X1 (using that Vi<^i is C a ), 
9k! £ C a (lR 3 ) for all a G (0, 1), and so, by standard elliptic regularity, 
u G C 2 ' a (R 3 ) for all a G (0,1). To recapitulate, for any x x G <Si(o;), 
the function x i— ► (/?i(s,xi) satisfies 

¥ >i(x,x 1 )='u il (a;) + i|a;| 2 (c~), ^ eC 2 '°(R 3 ), q6 (0,1). 

4 p| 

(3.42) 

Note that 

C = ZVi^(0,x0 = ZV^O). (3.43) 

We now apply the above to prove the existence of the pointwise limit 
of the integrand in ()3.36|) for fixed wG§ 2 and Xi G S x (uj). 
First, since ip = e~^^(pi, we have, for j — 2, . . . , N, 

[Vi(|V^| 2 )(ro;,x0] = u ■ [Vi( e - z|a:| |V^i| 2 )(ra;,xO] 

= -Ze- Zr \\/ m \ 2 {ru^ x ) 
+ e- Zr {u- [VidV^il 2 )^,^)]}. 

Because of ()3.2()j) in Lemma 13 .9| this has a limit as r j for fixed 
u G S 2 (recall that j = 2,...,N). The contribution to lim r ^ ^i( r ) 
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from this is 

N 



-2V / |V^i(0,xi)| 2 dxicL;, (3.44) 

~L i§2 J R 3N-3 
3 — z 

since terms proportional with to vanish upon integration. 

So we are left with considering uj ■ Vi(|Vi^| 2 ) (see ()3.36p ). To this 

end, use again ip = e~^ x ^\ to get 

w- V X .(|V^| 2 ) =u- V x {\ -f^ + e-f^v^l 2 ) (3.45) 
= u ■ V,(f^ 2 + e~ z ^\V m \ 2 - Z{oo ■ V^Oe-^Vi)- 

(We leave out the variables, (ra>,Xi)). We will study each of the three 
terms in (|3.45jl separately. 

For the first term in (|3.45|) . again using ip = e~^^ifi, we have 

" ' V.(f ^) = -fy>* + ^e'^(u ■ V m ), (3.46) 

which has a limit as r j for fixed w 6 S 2 and Xi G «Si (<*/)> since y?j is 
C 1,a (see (|3.25jl ). The contribution to lim r ^ t{(r) from this is 

-?/ / |^(0,x 1 )| 2 dx 1 dc;= -^Pi(O). (3.47) 

4 JS 2 JR3N-3 4 



As for the second term in ()3.45|) we have 



CO 



V.(e- z M|v lW | a ) = -Ze-^l^lVx^l 2 + e"^ ■ Vi|V lV i 



2\ 



(3.48) 



where the first term has a limit as r j for fixed w 6 S 2 and xi G iSi (a;) 
since ipi is C 1 '" (see (J3.25j0 . This contributes 



-Z / |Vi¥>i(0,xi)| 2 dxi<L; (3.49) 

JS 2 JrSN-3 

to lim r ^ 1{ (r). For the second term (in (|3.48jl ). using that <p\ = u + 
±r 2 (c • w) (see we get that 

|Vi^i| 2 = |VH 2 + -[(co- Viu)(c ■ x) + r(c • Viu)] 

3, ^9 If \ 

+ iq\ c ' x ) + 7g r ( c ' c )' 
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and so 

3 1 

w ■ Vi|Vi^i| 2 = 2(u;, (D 2 u)Viu) + ~r(c • u) 2 + -r(c • c) 

o o 

■ x) + (u) ■ Vi«)(c • cj) 
+ (c-V 1 u)+r(uj,(D 2 u)c}] . 

Here, -D 2 w is the Hessian matrix of u = with respect to x, and ( •, •) 
is the scalar product in M. 3 . Since ipi G C 1,a , and w G C 2 ' a , all terms 
have a limit as r j for fixed u G § 2 and xi G <Si(u;). We get 

lim [(e-^ ■ VilV^il 2 )^,^)] = 2u; ■ [(D 2 u^)(0)V xUil (0)] 

+ \ [(c • V A (0)) + (w ■ V^(0))(c ■ w)] . (3.50) 
When integrating, this contributes 

/ / |Vi^i(0,xi)| 2 dxidu; (3.51) 

JS 2 JR3JV-3 



2Z 
T 



to lim r |o (r); to see this, use ()3.43j) and Lemma lA. 21 in Appendix El 
and that J g2 k> = 0. 

For the third and last term in ()3.45|) . using that <£>i = u + \r 2 {c ■ to) 
(see (j3.42|0 and u ■ V x = d r , we get that 

u ■ V x .( - Z(u ■ V lVl )e- z \ x \y x ) = Z 2 (lu ■ Vi^e^l^ 

- Ze- z W<p x [ (u, {D 2 u)u) + \{c-u)]-Z{u- V x <pi) 2 e~ z ^ . 

Again, since (pi G C 1,a , and u G C 2,a , all terms have a limit as r | 
for fixed w6§ 2 and Xj G 5i(u;). 

The contribution from this to lim r jo t{(r) is, using Lemma lA.2| and 
Jg 2 u du = 0, 



3" 



3 
3" 



Vi^iOD,*!)! 2 ^!^ (3.52) 
A x Mx 1 (0)(^i(0,xi) dxtdui. 



§2 JR3N-3 



Here we used that Tr(D 2 f) = Af. 

This proves the existence of the pointwise limit of the integrand in 
(13~3oT) for fixed w6§ 2 and x x G 5i(w). Also, from (j3~l4l . (jSHj), (l3~4T^ 



THIRD DERIVATIVE AT THE NUCLEUS 
(JUlEJ), and ljTS2j) . 



21 



hint/ (r) 

r|0 



+ 




|V¥»i(0,xi)| 2 dxidw + — Pi(0) 
|Vi<^i(0, xl)| 2 dxx do; 

§2 JR3JV-3 




-— / / A x Mx 1 (0)v9i(0,x 1 )rfxic/u;. (3.53) 

•' Js 2 </R 3Ar - 3 



Note that, due to (pTT^I and ^ = e~T W^j, 

A^^O) = -A xl ^(0, Xl ) - + £- Viv-i,z-i(xi))V(0,xi) 
since Vi(0, Xi) = Vat_-i,.z-;l(xl) (see fll.3|) ). This implies that 



limtiYr) 

rj.0 



+ 




Z 3 

IV^i^xO^dXidw +— pi(0) 
'a 2 Jr^-* 4 J 

Z 2 

|Vi^i(0,x 1 )| 2 c/x 1 ^ + — pi(0) (3.54) 

/§ 2 jR3iV-3 4 





/§2 JR3JV-3 

Here we used that 



| V 4l V(0, Xl )| 2 + (Vat-i^-iCxO - £7)|V»(0, Xl ) 



dxi da; 



-/ ?/>(0, xi) (A Xl ■?/•>) (0, xi) dxi = / |V Xl ^(0,xi)| 2 d Xl . 

JR3JV-3 jR3iV-3 

(3.55) 

Before we provide a dominant in the sense of (j3.15j) to the integrand 
in ()3.36|) we now compute ii(0). Note that for r > 0, 



h(r) 




S 2 J R3JV— 3 



|Vi"0(ra;, xi)| 2 dxi dw 



(3.56) 



+ 52 / / |V^(^,xi)| 2 dxi^. 

~L J§2 JR3JV-3 



Use ip = e 2l x ly?i, then, for all w 6 § 2 and x x G <Si(a;) fixed, and 
j = 2, . . . , N, Lemma 1331 gives that 

\V j iP(ru,± l )\ 2 ^e- Zr \V j( p l (ru,± l )\ 2 ^ |V^i(0,xi)| 2 . (3.57) 
In particular, this proves the existence of 



N 



N 



|V Xl ^(0, Xl )| 2 = |V^(0, Xl )| 2 = |V^i(0, Xl )| 2 . (3.58) 
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Similarly, 

Z 2 

|Vi^(rxu,X!)| 2 = — |^(rw,Xi)| 2 + e- Zr |Vi^i(rw,Xi)| 2 

- Ze~^ r i/;(ruj, ii) (u ■ Vi^i(rcj, Xl )) (3.59) 
^ ^|V(0, xi)| 2 + | V a¥ >i(0, xx)| 2 - ^(0, *0 (a; ■ V^O, x 2 )) . 

Using again L 2 udu = 0, it follows from (J3.56j) - (J3.59j) . and Lebesgue's 
Theorem of Dominated Convergence that 

limti(r) = ^-jBi(O) + / / |Vi^i(0,xi)| 2 dxidu; 

+ 11 IV^O,^ 2 d±iduj. (3.60) 

This proves ([3.9|) . Combining this with (J3.54|) (using (J3.58jl ) proves 

Dominant. We turn to finding a dominant to the integrand in ([3.36)1 . 
We shall apply results from j^j, recalled in Appendix IE1 

From the a priori estimate (|B.8|) in Theorem IB.3I and (J3.36)) follows 
that, for almost all (x, xi) = x £ M. 3N (choose, e.g., R — 1, R' — 2) 

u ■ V 1 (\ViP\ 2 ){x,± 1 )\ < C |V^(x, Xl )| U\\ L ^(B 3N (( X M),2)) (3-61) 

+ |fVVV2 Xl ' fc f d ^ \j> ° 2Fcut )(x xi) 

with F cut = F 2iCUt + F 3)Cat as in Definition IB. 21 

Now, using the exponential decay of ip ([1.8)1 (assumed), and of 
(|1.9j) . we get that there exist constants C, 7 > such that 

|V^(a;,xi)| 1 1 -0 1 j z,oa (S3jv((ar,* x ) ,2)) < Ce~ 7|xi1 for almost all x £ M 3 . 

(3.62) 

This provides a dominant (in the sense of ([3.15)1 ) for the first term in 

(jmji. 

We need to find a dominant for the second term in (|3.61jl . First 
recall that F cut = F 2)CUt + F 3 cut . With F 2 as in ()B.1|) we have F 2iCU t = 
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F 2 + (F 2iCUt -F 2 ), with 

N Z 



(F 2 , cut - F 2 )(x) = E " 2 (x(kil) " 1) W (3-63) 

i=l 

+ E tWI^-^I)- 1 )!^-^!- 



4 

l<i<j<N 

Note that <9 2 [(x(M) — is bounded in M 3 for all second derivatives 
d 2 , due to the definition (jB.3|) of x- Using the exponential decay of ip 
(II. 8|) . and (|1.9|) . we therefore get a dominant for the term 

V V V 2^ f ^ V d 2 ( F 2,c U t -i^h ( x x x ) 
\xi\\dx l>m ) dx ltk dx e , m / ' 1 ' 

f=l fc=l m=l ii> 11 

A tedious, but straightforward computation gives that 

yyy 2 ^.(d±\ i , (3 . 64) 

£=1 k=l m=X 
N 



5*E 



i=2 



1 xi 



V \x\\ \x\ — Xir / L v 



We first remark that, again using exponential decay of i/j and V^>, we 
get the estimate 

V / . \X\ XA \X\ / 
1=2 1 



N -c\xi\ N 

- c ^[<tist(L(u),x i ) ( n e 

i=2 V V h 11 3=2 j& 

where dist(L(u;), y) is the distance from y E M 3 to the line L(co) spanned 
by uo E E> 2 . Note that for fixed uo E § 2 , the function e~ c \ yl / dist(L(cu) , y) 
is integrable in R 3 , and its integral is independent of u E S> 2 . Therefore 
the right side of (|3.65|) is integrable in M 3JV_3 for fixed uj E S> 2 , and is 
integrable in IR 37V ~ 3 x § 2 (all uniformly for r E M). 

The argument is similar for the second term in ([3.64)1 . 

We are left with considering the terms in (j3.61|) with i*3 )CU t. Let 

f 3 (x,y) = C Z(x-y) \n(\x\ 2 + \y\ 2 ) 
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(so that i*3(x) = J2i<j h{ x ii x j)i see (IB ■ 2 j) ) . For all second derivatives 
d 2 we easily get, due to the definition (|B.3|) of x, 

d 2 [x{\x\)x{\y\)h(x, y)] = X(\x\)x(\y\)d 2 f3(x, y) + g 3 {x, y) 

= C Z X (\x\) X (\y\) ln(|x| 2 + \y\ 2 )d 2 (x ■ y) + g 3 (x, y) , 

with (73,^3 bounded on R 6 . Therefore, defining, for all second deriva- 
tives d 2 and with x as above, 

^F 3 ,cut(x) := C Z J2 X(l^l)x(kil)ln(k| 2 + |y| 2 )a 2 (x-y), (3.66) 

l<i<j<N 

we get a dominant for the term 



N 3 3 



Z / Z , / Z , / r -, 



■0 



<9 2 F- 



3, cut 



<9 2 F 



3, cut 



dxiudx. 



£,m 



dxiudx 



£,m 



IX, *i 



l=\ k=l m=l 

using the exponential decay of ip and V^. We find that 

N 3 3 



£=1 k=l m=l 



xi, k ( di) 



\xi\ \dx 



d 2 F 



3, cut 



dx\ kdx 



N 



2C Z^J2x(\xi\)x(\x l \) H\ Xl \ 2 + 



x ; 



i=2 



X\ 
\x\\ 



(3.67) 



Note that, by the definition of Xi 

X {\x\)x{\y\) | ln(|x| 2 + \y\ 2 )\ < X {\y\) [| ln(|y| 2 )| + 3] 
and so, again by the exponential decay of ip and Vip, 

N 



8=2 



2C ZijJ2\x(\xi\)x(\x l \)H\x 1 \ 2 + \x l 



N 



Xl 

\ X l\ 
N 



[ruj,x 1/ 



(3.68) 



<G(E^N)0 ln (N 2 )l+3])(n e " cW 

i=2 j=2 

for all ijj G S, r G K and (almost) all Xi G «Si(o?). This provides a 
dominant for the term 

N 3 3 



=1 k=l m=l 



<9 2 F 



3, cut 



dxikdx. 



Lm 



(X,Xi 



and we have therefore provided a dominant, in the sense of of (|3.15|) . 
for the integrand in (j3.36j) . 

This finishes the proof of Proposition 13.31 □ 
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Appendix A. Two useful lemmas 

The following lemma is Lemma 2.9 in we include it, without 
proof, for the convenience of the reader. (The proof is simple, and can 
be found in [S]). 

Lemma A.l. Let G : U — > lR n for U C R"+ m a neighbourhood of 
a point (0,?/o) G R™ x M m . Assume G(0,y) = for all y such that 
(0,?/) e U. Let 



Then, for a G (0, 1], 



x = 0. 



G e C°' a (U; R n ) f e C°' a (U). (A.l) 
Furthermore, \\f\\c<*(u) < 2\\G\\ c <*(u). 

The following lemma is used to evaluate certain integrals. 

Lemma A.2. Let a, b e R 3 and let A <G M 3x3 (C). 
Then 

4ir If 
(uj ■ a){uj ■ b) duj = —{a ■ b) = -(a ■ b) I du , (A.2) 

■2 3 3 ./;;■ 

/ to ■ {Aw) du = —Ti(A) = -Tr(A) [ dcu . (A.3) 
7s 2 3 3 J S 2 



Proof. Both (|A.2|) and ()A.3|) follow from the identity 

/ UiUJj duj = — 6ij , (A.4) 
Js 2 3 

which we now prove. 

For i ^ j, the integrand UiUj is odd as a function of Uj, which implies 
flA.4|) in that case. For i = j we calculate, using rotational symmetry, 

/ ufdu) = ~ I (uj\ + ul + ul) du = 
Js 2 <J Js 2 <J 

This finishes the proof of (|A.4|) and by consequence of Lemma [A. 21 □ 



Appendix B. Regularity of the eigenfunction t/> 
The following two theorems were proved in jHj. 
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Theorem B.l. (jHl Theorem 1.1 for atoms]) Suppose ip is a solution 
to Hip = Eijj in Q C R 3N where H is given by (JT2J). Let T = e F2+Fz 
with 

N Z 1 

*2(x) =^ -— 5Z ll^-^'l' ( R1 ) 
i=l l<i<j<N 

F 3 (x)=C ^ lnd^ + l^f), (B.2) 

where Co = 4=^. T/ien ^ = JF0 3 witt 3 e C 1 ' 1 ^). 
Definition B.2. Let \ e CfW, < X < !> with 





for x 


I" 


for x 



We define 

+ F 3)CUt , (B.4) 



-Tcut — -T2,cut T -T3,cut) 

where 

? Z 1 

^2 )C ut(x) = 53 - 2"x(kii) kti + ^ 5Z x (' x * ' Xi ( R5 ) 

i=l " l<i<j<N 

F 3 , cut (x) = C 53 Zxd^Dxdxj-D^.Xj-) lndx^ + lx.f), (B.6) 

l<i<j<N 

and where C is the constant from (jB.2|) . We also introduce 3iCUt by 

V> = e Fcut 3 ,cut. (B.7) 

Theorem B.3. ([BJ Theorem 1.5 for atoms]) Suppose ip is a solution 
to Hip = Eip in M. 3N . Then for all < R < R' there exists a constant 
C(R,R'), not depending on ip nor x G M. 3N , such that for any second 
order derivative, 

&= * d \ , i,j = l,2,...,N, k,£ = 1,2,3, 
OXi tk dx jt e 

the following a priori estimate holds: 

\\d 2 iP-ipd 2 F cut \\ Lac(B:iN ^ hR)) < C(R,R!)U\\ La o (B31f(:KOtRf)) . (B.8) 
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